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We proposed a new formulation for heavy quark effective tlieory (HQET), wfiose Lagrangian is 
fiermitian and has a manifest reparametrization invariance. As an application, we calculated the 
semileptonic and nonleptonic inclusive heavy hadron decay rates up to second order mass corrections 
and found that there are no kinetic energy correction terms. 
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^ ■ I. INTRODUCTION 

^ ■ o 

PsJ ' The heavy quark effective theory (HQET) is an useful tool in studies of heavy quark systems. In infinite heavy 
hj^. quark mass limit, there exists a heavy flavour-spin symmetry of QCD. To derive a HQET theory from QCD, one 
' needs to integrate out the high frequency modes of the heavy quark field to retain the low frequency modes. The 
, high frequency modes with scales larger than two times heavy quark mass can be integrated out by employing the 

• equation of motion (EOM) method ||^ or the functional integration (FI) method Q . The Lagrangians of these two 
T-H ' HQET theories contain non-herniitian terms which would lead to imaginary mass corrections [Q. Nevertheless, this 
04 ■ problem did not receive too much attentions in literature. We will show in Section II that the non-hermitian terms 

can be regularized by further integrating out the high frequency modes with scales larger than heavy quark mass. 
. The heavy quark momentum can be factorized into the heavy quark velocity part and the residual momentum part. 
Reparameterizing the heavy quark velocity and the residual momentum would not change the effective theory. It 
implies that the coefficients of mass correction terms in HQET Lagrangian can be fixed by reparameterization 
It was noted that the EOM and FI theories have no manifest invariance under reparameterization 1^0]. Another 
lO HQET theory was derived by employing the Foldy-Wouthuysen (FW) transformation It is unclear that whether 
the theory is invariant under reparameterization. 
, ^ The non-hermitian terms and lacking a manifest reparameterization invariance should be avoided for a self-consistent 

HQET theory. We will show that a new formulation for HQET could simultaneously resolve these flaws. The remainder 
of this paper is organized as follows. Section II will explain the cause that leads to the non-hermitian terms, and 
Oh, Section HI will present the derivation of a hermitian Lagrangian from QCD. In Section IV, we will show that the 

• hermitian Lagrangian is manifestly invariant under reparameterization. Section V devotes to evaluation of mass 
^ ' corrections for inclusive B decays. An Appendix enumerates the mass correction terms up to 0{1/Mq). 

^ • 



> 
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II. THE NON-HERMITICITIAN TERMS 

The happens of non-hermitian terms in EOM and FI theories could be better understood by investigating the 
non-relativistic reduction of the Hamiltonian of an electron interacting with electromagnetic fields. For simplicity, 
we consider the case that the EM fields are static. The equation of motion for an electron under static Coulomb 
potentials takes the form 

= [a ■ TT + (3m + eV]'ilj , (1) 

ot 

where V represents the Coulomb potential, m denotes the electron mass, is the electron wave function, {d)i ~ ^'^^^ 
with i = 1, • • • , 3, /3 = and tt = — iV. In the nonrelativistic limit E ^ m + jy^/2'm, it is convenient to use the large 
(j) and small % components of ip defined as 



to recast (0) into two coupled equations 



^^il) (2) 

A 



i— — a ■ TTY + eVcj) + mcj) , 
ot 

i— = a-Tf(f) + eVx-mx, (3) 
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where we have employed the representation 

As times evolute, the potential contributions are smeared out by the large value of electron mass, m. To avoid this, 
one may introduce slowly varying functions of times $ and X defined as 

$ = e'™V , X = e™*x, (5) 

whose equations are easily derived 

i—- = a ■ ttX + eV^ , 
at 

dX 

i— =a-7r<i> + eVX -2mX. (6) 
Because that eV <C 2m, we are legal to approximate X into 

X = — -a ■ 7?$ 

2m + vr" 



2to 4m-^ 

and substitute the expanded X into the first equation of (0) to obtain 



+ •••]*, (7) 



• — = eVS} + '— ^ 

dt 2m 4m^ 



where 

[(? • TTV]a ■ 7? = Try • 7? + i(7 • (ttF X tt) . (9) 

After rewriting (H) into Hamiltonian 

H = ^^eV^ + - ^I'-^jfa • nV]a ■ tt + • vf)^}^ , (10) 

2m 4m2 J v / j > v y 

we note that the Darwin term in the above Hamiltonian 

Or, = -%$^7fF • 7?$ (11) 
4m^ 



is nonhermitian 



0^^ = ^{^^^ ■W)<^^Od ■ (12) 



One can add up 0|) and Od , and divide the sum by 2 and perform integration by parts by ignoring the surface terms 
to derive a hermitian Darwin term IBtl 



og = £^$nw'^]'i> • (13) 

Or equivalently, one can make use of the renormalized $ with expression 

$jvfl= l + ^ + --- * 14 
8m^ 

to obtain the regularized Darwin term Og. The regularized Darwin term Og is in fact the second spatial variations 
of V due to the jittery motions of electron with Compton wavelength &r \/m 

< V{r + dr) >«< V(r) > < {nfV > , (15) 

om^ 
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where the bracket means integration with electron wavefunctions and the first order variational term vanishes due to 
the assumption that the wave functions are spherically symmetric. 

The above example exhibits the cause of the nonhermitian terms. Similarly, we can show the non-hermitian terms 
in the HQET theories derived by EOM or FI. The equation of motion for the heavy quark field "0 is 

{ip - Mq)^ = , (16) 

where Mq denotes the heavy quark mass and ip is the covariant derivative ip = ip — g/lT" . At energies much 
below than Mq scale, tp is not an appropriate variable. One invokes field redefinition Q{x) = exjp {iMgv ■ x)4'{x) 
to remove the large phase factor Mqv from the wave function. The v variable represnts the heavy quark velocity. 
Rewriting ([l6| ) in terms of Q{x) yields 

{tp - 2AfQ^i^)Q - . (17) 

By imposing condition = 1, one can separate Q into large h and small H components 

1 + 1^ 1 — /' 

Q=^Q + ^Q = h + H . (18) 
Substituting (|l|) into (|l^) and multiplying (1 - ^)/2 from the left of (^ yields 

g= ... \ -n i^P±)h (19) 
with = V ■ D and p± = p - fD\\ . Using (|l|) and (|l|) leads to 

Q = [l+ ^,, i^Pi-)]h. (20) 

2Mq + iU\\ 

Substituting (|^) into (|l|) and expanding it up to 0(1/A/q), one then arrives at 

-^[^'^o.pvxG'^^Di + iD\\a^pG"P - iD\\Dl + v^G^^ Dj,] + 0{^)]h , (21) 

where ^''7'^ = + ia^'^ and [il?^, iD'^] = —iG^'^ have been used. Note that the Darwin term (the last term in the 
second line of (pTj)) is nonhermitian. Following the route in previous QED example, we may regularize the Darwin 
term by renormalized large components h' 

h'^{l + ^^Pl + (22) 
The equation of motion for h' thereby takes the form 

iD^h' ^ \ —\-D] +-a-G] 

" ^ 2Mq ^ ^ 2 ^ 

+ ^[^'^o.fiV,{Dl, G^'^} + Va[Dj,G"0]] + 0{^)]h' . (23) 

The Darwin term in the second line of (^) correspond to the relativistic effects of Zitterbeweguen from the jittery 
motions of the heavy quark with Compton wavelength Xq « 1/Mq. This implies that the large components h still 
contains high frequency modes whose scales are larger than Mq. These high frequency modes should be integrated 
out to have a low energy effective theory. It concludes that the mass corrections receive two kinds of contributions: 
the first kind of contributions comes from the integration out of heavy antiquark degrees of freedom, and the second 
kind of contributions arises from the integration out of the high frequency modes of the heavy quark degrees of 
freedom. Only both kinds of contributions having been carefully taking into account can result in a hermitian th eory . 
The integration out of the high frequency modes is equivalent to the renormalization for the heavy quark field (p3). 
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This means that the renormalized field h! contains only low frequency modes with scales not larger than Mq and is 
shwon responsible for low energy physics. A systematic method, which can derive a hermitian Lagrangian as well as 
an appropriate effective field, is very intriguing in theory and phenomenology. To develop this method is the main 
purpose of this paper. 

To reveal the eligibility of the unrenormalized large components h, we support two examples. The first example is 
the spin sum of h in free theory. By definition (^0|), the spin sum takes expression 



Y^h{X)h{\) = i±i^Q(A)Q(A)i±^ , (24) 



A A 

where A denotes the spin of the summed spinors. The spin sum over Q is equal to 



^Q(A)g(A) = i±^ + -^, (25) 



A - --'Q 



2 2Mq 

where k means the residual momentum whose magnitude is much smaller than Mq. Substituting (|25|) into ( ^4|) yields 

rf) (26) 



y > y I 2 ^ 4M/ 

A Q 



It is noted that the spin sum of the HQET effective spinor hy is equal to 



Y,K{\)K{X) = . (27) 

A 

This example shows that the h propagator differs from the /i„ propagator 

- (-) 

by a factor (1 — -^^), which is just twice the reversal renormalization factor in (^). We take the free heavy quark 
spinor as the second example . Let uq denote a free full heavy quark spinor whose energy is Eq and mass Mq and 
spatial momentum k. uq can be expressed in terms of its rest frame spinor as 



= I ; IT ... I (29) 

yf2MQ(MQ+EQ) 

where 

^'''-{l) ' '^'''-{\) (30) 
denote the rest frame spinors. In the static approximation, we can expand Eq 

Eq^^M^ + P^Mq-^ , (31) 
where P = f^nd fc^ = (0, k) have been used. Under this approximation, the full spinor uq becomes 



From (|32D, we can identify the large components h and small components H of uq 



"^^f-%U^^^"V- '''' 



2Mq ~ ft± 

(33) 
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Spinors (p and uq are well normalized (j)^(f> = uquq = 1, while the large components h has an incorrect normalization 
hh = 1— {^±/2MqY as pointed out in Finally, wc emphasize that from the equations of motion of Q (|l9|) we can 
directly derive the relation between h and H as 

and on-shcU condition for Q 

III. CONSTRUCTION OF THE EFFECTIVE THEORY 

A. Derivation of The Effective Field 

One may match QCD to HQET at the heavy quark mass scale by requiring that the IPI Green's functions of the 
two theories describe the same physics. The simplest way to achieve this is by setting the external quark momenta 
on shell |^]. The LSZ reduction formula for a heavy quark fermion in momentum space is 

SiPQ, ...) = ^lZ^(FQ/MQ)i^— ^ . . . / dxe^^^- < 0\Tmx) . . .]|0 > \p. 



Q{v + k/MQ){^ -K)-- - I dxe'^-- < 0\T[Q{x) . . .]|0 > U=-ky2MQ (36) 



, where Q{x) = exp [iMqv ■ x)'ip{x) and ipi^) denotes the heavy quark field. To proceed, we employ the projection 
operator method to derive the matching between the effective spinors Q and hy. Q carries momentum Pq = Mgv + k, 
while hy has momentum k with respective to a constant moving frame with velocity v. Both spinors are equivalent 
variables for low energy physics. However, in the limit Mq ^ oo or at the energy scale much below Mq, hy is more 
appropriate than Q. 

In this method we want to develop, we can specify Q by a positive energy projection operator 

which selects spinor Q carrying momentum Mqv + k 

A+Q ^ Q . (38) 

Eq.(^8|) is equivalent to the equation of motion [ip — Mq{1 — 'fj)]Q{x) = 0. Being a projection operator, A+ obeys 
the identity 

(A+)2 = A+ (39) 

, which implies on-shell condition 

[fci, + = (40) 

with fc|| = V ■ k. The inverse operator of A"*" is the negative energy projection operator 

A-^(l^-^ (41) 

2 2Mq ^ ' 

defined by identity 

A+ + A" = 1 . (42) 
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In order to derive hy which respects the physics in the Umit Mq — > oo, we define projectors 

= ^ ^ ^ lim A± . (43) 

Note that operators A+ (A~) arc the infinite mass hmit of the energy projection operators A+(A~). A^ satisfy 
identities 



By 1 = A+ + Ay , we recast Q to be 



From (|38D and (45) we arrive at 



and 



{Atf^At- (44) 
Q = A+Q + A-Q = h + H. (45) 



^ ]h . (47) 



In hterature, people always stop at this point and identify h as hy. However, as point out in last Section, h is not 
identical to hy. It is natural to assume that hy is the limit Mq — > oo of /i and two spinors are proportional. The first 
assumption comes from hy = limMg— >ooC? , while the second is based on the fact that both h and hy are projected 
out by A+. In this way, we argue that h — [1 + with ansatz to = u and '^uj — uj. To derive w, we note identities 

A+-^QQ, (48) 
A+ = ^ hyhy , (49) 

where summations over spin indices of the spinors are implied. Eqs. (|4^) and (^) hold if and only if QQ = hyhy = 1. 
Eq. ( ^9|) c omes directly from definition, the limit Mq — > oo of (^8|). Substituting (^ ) and h — [\+uj\hy into (|4^ ) and 
using ([49|) for hy, we get the equation for ui 

(^)2 + 2c^ + i±^ = (l-i)(i±^ + ^)(l-^) . (50) 

With the help of on-shell condition (pO|), ( pO| ) is recast as 

("'' + ^" + <2]fe'"<^' = »- 

The above equation is easily solved 

uJ = -l + ^/lTT (52) 

with T = -{2M^?{^)- We then obtain 



The relation between h and /i^, is consistent with that between h and h' found in last Section. 
Combining (47) and (p3|), we get 



Q = J\^^^j^^^Ayhy=A{w = v + k/MQ,v)hy. (54) 
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Note that ( p4| ) is just the Lorentz transformation between two spinors with relative velocity k/MQ. The transformation 
operator K{w = v -\- k/Mq, v) is identical to the Lorentz boost in the spinor representation Q 

KM^^m= ,55) 

^2(1 + v ■ w) 

In presence of interactions, the Lorentz boot interpretation of (jsj) is invalid. The reverse transformation from into 
Q can be derived in a similar way. The result is 



_ / l-|^/(2AfQ) 



Transforming (p4) into coordinate space by replacements f, ip^ one obtains 



Field Q{x) is consistent with the field derived by Luke and Manohar [|| 

^^{x) = K{v + iD/MQ,v)K{x) . (58) 
By employing (|5^ ) and (^7|), we arrive at the matching between Q and /ly at scale equal to Mq 



S{k,...)^^h,{v){-l--l-^)... I dxe^''-^0\T[K{x)..mU^-ky2M^ , (59) 



which is different from (41) in ||10|] by a factor y Z(fc) = w(l — ft'^/{4:MQ). This is because the effective field /i„ 



employed in po| is the unrenormalized large components h^'-' ~ h = y Z{k)hy. 

By matching QCD and HQET at 2PI and quark- gluon-quark interaction Green functions, we can derive the HQET 
Lagrangian 

L = i^{ip - Mq)^ 
= Q{tp^2MQA-)Q 



— + 1+^P/(2Mq) 1+^P/(2Mq) 

= ^^^^ V i-zp/(2Mq)(^^ - ) V r^jprn^) ^ ^ ■ ^ ^ 

Note that the HQET Lagrangian is hermitian. This can be inspected by explicitely performing hermitian conjugation 
and integration by parts. 

IV. THE VELOCITY REPARAMETERIZATION TRANSFORMATION 

A. Field Transformation 

The heavy quark momentum Pq makes no differences for parameterizations Pq ~ Mqv + /c or Pq = AIqv' + k' . 
It was found that the HQET Lagrangian should be invariant under reparameterizations for the velocity w — > and 
for the residual momentum fc — *■ A;'. We show the following theorem: If v and v' relate each other as, v' = v + 5v with 
(ti')2 — ■i;2 — 1 and V ■ 5v + ((5u)^/2 = 0, then h^i relates to hy as 



and 
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hA^) = ^^''"'''■yT^f7^A^^'"(^) ■ (62) 

Furthermore, if MqSv ~ k — k' , Q{Pq = Mqv' + k') equals to Q{Pq = Mqv + k) and 

Q{Pq = Mqv' + k'){x) = e^*^«*'^-"Q(PQ = Mqv + k){x) . (63) 

The proof of this theorem is straightforward. Since v' — v + 5v and (v')^ — = 1, velocities v' and v have 
corresponding energy projectors (1 + ^')/2 and (1 + ji) /2, which will project effective fields hyi and hy, respectively. 
By replacing 5v with k/Mg in the transformation (p3) , we thus derive the transformation from hyi to hy. The proof 
of Q{Pq = Mqv' + k') identical to Q{Pq — Mqv + k) is also trivial by noting that 



= + = Ji±M|MAj,... ,65, 



and MqSv =^ k ~ k' . It results in the identity 



' 1 + |^7(2Mq) / l + J^/2 
1-^'/{2Mq)V 1-5^/2 " ^ 



1+1/(2Mq) 

^Q{Pq^ Mqv + k) . (67) 

The proof is completed. The importance of this transformation is the association of successive transformations. If we 
denote the transformation from hy into hvi=v+sv by hyi = L{v,v')hy, then we can have L{v,v") = L{v' ,v")L{v,v'). 
We show this explicitely below. The successive transformations v ^ v' = v + Svi followed by v' v" = w' + 5v2 = 
V + Svi + dv2 would have the transformations in the effective field as 



h^^hy, = J^-^^^h: (68) 



and 



' l + (3^1+ 6-^2) /21 + f' + 
l-(<5^i+<5^2)/2 2 ^ 



(70) 



B. Reparameterization Invariance 

We now show that the reparameterization invariance is trivial and manifest for the Lagrangian (|6^). By previous 
theorem, it is also straightforward to prove that the effective Lagrangian in terms of Q is invariant under transforma- 
tions V v' = V + Sv, MqSv k - k' and Q+{Mqv' + k'){x) = exp {IMqSv ■ x)Q+{Mqv + k){x) 



8 



L+ = Q+{Pq = Mqv + k){ip- 2MqA-)Q+{Pq = Mqv + k) 
= Q+{Pq = Mqv' + k'){ip - 2MQk~,)Q+{PQ = AIqv' + k') 

It is also trivial to prove the above invariance for effective Lagrangian in terms of 

L+ =. Q+{Pq = Mqv' + k'){Up - 2MqA^,)Q+{Pq = Mqv' + k') 



(71) 



7+.+ l + ^P/{■2MQ)^^ 



1 - ipl(2MQ) 



{iP^2MqK-,). 



' 1 + iP/{2Mq) 
1 - iP/{2Mq) 

{ip - 2MqA 



-^''^^P-S^/2' \^l-^p/{2MQf^ ^^^''^^'>^l-^p/{2MQf P- 



-6^/2 ^ " 



= Q+{Pq = Mqv + k){ip- 2MqA-)Q+{Pq - Mqv + k) 
, where we has used 



/ 1 + ^P/(2Mq) / 1 + W2 _ I i + iP/{2Mq) 



1 - iP/{2Mq) 



1 - <5/V2 V 1 - ^PH'^Mq) 



(72) 



(73) 



C. Zero recoil limit 



The matrix element of {H.t,\QTQ\H^) appears much simplicr in the zero recoil limit, v' = v. With (|57| ) and (|59|), 
we can derive the following identities 



{H,\QQ\H,) 
{H,\Ql^Q\H,) 



{Hy\hvhv\Hy) , 



(i/.|/i.(w^ + — )/i,„|77,) , 
^P 



Mn 



{H,\Qa^''Q\H,) 
Notes that the pseudoscalar currents vanish. 



Mq 



(74) 
(75) 

(76) 

(77) 



V. APPLICATIONS 



We apply the HQET Lagrangian (|60|) to calculate the rates for inclusive semileptonic and nonleptonic heavy hadron 
decays up to second order mass corrections. According to operator product expansion (OPE) [p^|-p^, the decay rate 
of a heavy hadron Hb containing a b quark are expanded in terms of local operators with increasing dimensions 



{Hb\cibb + ccb-crGb 



\Hh 



(78) 



where Ci and cg are short distance coefficients. The momentum carried by Hb is chosen as Ph^ = Mh^v with Mh^ 
the Hb mass and v its velocity. \Hb) means the eigenstate of the full HQET Lagrangian (^0|) with normalization 
{Hb{P')\Hb{P)) = v^ {2ttY 5'^ {v' — v). The next step is to expand the matrix elements of the local operators into 
inverse powers of the &-quark mass. Because of ( [74[ ) the matrix element {Hb\bb\Hb) is transformed into {Hb\hyhy\Hb) ■ 
There is no mass corrections in {Hb\bb\Hb) . This differs from the conventional result 



{HbmHb} = i- '''^^'''\-f-^'''K oi^: 



2ml 



(79) 
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where is^^Hb) and /i^(iJb) parameterize the matrix elements of the kinetic-energy and the chromo-magnetic operators, 
respectively. The matrix element {Hi,\b^aGb\Hi,) has coefficient of 0(l/m^) and its leading mass corrections are of 
0{l/ml) will be ignored. With these considerations, we derive the semileptonic decay rate 

which is free from /i^ uncertainty in contrast with the usual formula 

- - - '"'-^W^.f (.o(^)(i - 4r^) - 2.,(!^)4 + ...). (81) 



1927r'^ 1^ TOfc 2m^ nib 

The difference is due to different accounts for the time dilation in the Fermi motion of the heavy quark in the heavy 
hadron rest frame. Recall that the weak decay widths of muon 

r(A^--^) = ^-o(4)- (82) 
lOzTT'^ mf, 

One can see that our result is more reasonable because the leading term {Hb\bb\Hb) is normalized to unity in our 
HQET theory. We can derive the nonleptonic decay rate 

= ^\V,b\U.uNc {w,zo{^) - 2w,{z,{^) + Z2(^))4 • • -I ' (83) 
1927r-' I nib mb mf^ ' 



which is different from the conventional result 
The Wilson coefficients are defined as 

2,2, ClC2 C1C2 , . 

with ci = (c+ + C2)/2 and C2 = (c_|_ — c_)/2 and the phase space factors are 

zo{x) = l-8x + 8x^ - x'^ -Ux^logx , zi{x) ^ {1 - x)^ , Z2{x) ^ {1 - xf . (86) 



VI. CONCLUSION 

We have regularized the non-hcrmitian terms in the EOM type HQET theories Q-jl] up to 0{1/Mq). We found that 
the large components of the heavy quark field are not appropriate variables beyond leading order mass corrections. 
Only the renormalized large components, whose high frequency modes have been integrated out, can be a relevant 
effective field for low energies. In terms of the renormalized large components, the HQET Lagrangian (|6^) is hermitian 
to all orders in 1/Mq and contains manifest reparameterization invariance. 

Acknowledgments: We are grateful for helpful discussions to Guey-Lin Lin and Hsiang-nan Li. This research 
was partly supported in part by the National Science Council of R.O.C. under the Grant No. NSC89-2811-M-009-0024. 

APPENDIX A: MASS EXPANSION LAGRANGIAN 

We discuss the mass expansion of the HQET Lagrangian derived in (|60|). The HQET Lagrangian L is expanded 
into mass correction terms as 

00 r (ri) 

l-Etttt-^ (A1) 



,^„(2Mcj)" 

where the first leading terms L^"\n — 0, 1, 2, 3 are enumerated as follows 
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hv 



-2iDl + -{va[D0,G''^]+iac,0Vx{D^,G^"}) 



(A2) 
(A3) 

(A4) 



(A5) 
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